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Abstract 

Motivated by experimental observations of exotic free surface standing wave patterns in 
the two-frequency Faraday experiment, we investigate the role of normal form symmetries 
in the associated pattern selection problem. With forcing frequency components in ratio 
m/n, where m and n are co-prime integers that are not both odd, there is the possibility 
that both harmonic waves and subharmonic waves may lose stability simultaneously, each 
with a different wavenumber. We focus on this situation and compare the case where 
the harmonic waves have a longer wavelength than the subharmonic waves with the case 
where the harmonic waves have a shorter wavelength. We show that in the former case a 
normal form transformation can be used to remove all quadratic terms from the amplitude 
equations governing the relevant resonant triad interactions. Thus the role of resonant 
triads in the pattern selection problem is greatly diminished in this situation. We verify our 
general bifurcation theoretic results within the example of one- dimensional surface wave 
solutions of the Zhang- Vinals model [jjj of the two-frequency Faraday problem. In one- 
dimension, a 1 : 2 spatial resonance takes the place of a resonant triad in our investigation. 
We find that when the bifurcating modes are in this spatial resonance, it dramatically 
effects the bifurcation to subharmonic waves in the case that the forcing frequencies are in 
ratio 1/2; this is consistent with the results in Q. In sharp contrast, we find that when the 
forcing frequencies are in ratio 2/3, the bifurcation to (sub)harmonic waves is insensitive 
to the presence of another spatially-resonant bifurcating mode. This is consistent with the 
results of our general analysis. 

PACS numbers: 47.54.+r, 47.20.Ky, 47.35. +i 
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1 Introduction. 



Exotic free surface standing wave patterns, parametrically excited by two-frequency forcing, 
have attracted considerable attention in recent years, both experimentally [Q, ||, ^, ^) and the- 
oretically ^, @, ||, In this system, the surface waves are excited by subjecting the fluid 
layer to a time-periodic vertical acceleration with two (rationally-related) frequency compo- 
nents. This corresponds to a modification of a classic hydrodynamic problem, which dates back 
to observations of Faraday [[To), in which the surface waves are parametrically excited by a 
purely sinusoidal vertical acceleration of the fluid container [fi"l| . Triangular patterns Q , quasi- 
patterns [s), and superlattice patterns || are among the more exotic states that have been 
observed in laboratory experiments employing two-frequency forcing. Here, "quasi-patterns" 
are patterns with twelve- fold symmetry and are the hydrodynamic analogues of two-dimensional 
quasi-crystals. These patterns are not spatially-periodic, but do have long-range orientational 
order; their spatial Fourier transform exhibits twelve prominent equally-spaced peaks that lie 
on a circle. Although such patterns have been observed in single-frequency Faraday experi- 
ments HH 0, they are more readily observed with two-frequency forcing and examples 
have been seen with excitation frequencies in ratio 4/5, 4/7, 6/7, 8/9 ||. In contrast, "superlat- 
tice" patterns are spatially-periodic with structure on two disparate length-scales. They have 
been observed in the two-frequency Faraday system with a forcing frequency ratio of 6/7 ||. 
The observations of both quasipatterns and superlattice patterns have indicated that the pat- 
terns are synchronous with the forcing frequency. The Faraday wave problem is an attractive 
experimental system for studying the fundamental mechanisms behind the occurrence of such 
patterns because of the fast time scales involved and the number of tunable control parameters 

One motivation behind some of the early experiments on two-frequency parametric exci- 
tation of Faraday waves was to destroy the Zi symmetry associated with subharmonic waves 

||, [|. Specifically, in the classic Faraday experiment with single- frequency forcing cos(wi), 
the onset surface waves respond subharmonically with frequency u>/2 (Harmonic response 
is also possible for very shallow layers of viscous fluid [|l8| |l^] . ) In the situation of subharmonic 
response there is a discrete time translation symmetry t — > t + that is broken by the state 
of the system. This symmetry-breaking manifests itself in the governing pattern amplitude 
equations by suppressing all even terms. Whether or not this extra Zi symmetry is present 
has a particularly profound effect on the formation of hexagonal patterns |2G|; for instance, in 
many systems where it is absent, hexagonal patterns are preferred at onset |21j| . Edwards and 
Fauve H noted that by introducing the second perturbing frequency component to the periodic 
forcing, they could destroy the discrete time translation symmetry of the system. 

A second feature of using two frequency forcing is that it is possible to obtain a neutral 
stability curve with minima at two distinct wavenumbers, where the ratio of these two crit- 
ical wavenumbers can be adjusted by varying the two frequency components of the forcing. 
Typically, one minimum corresponds to standing waves that are subharmonic with respect to 
the forcing frequency, while the other minimum corresponds to synchronous waves |Q. It was 
proposed by Edwards and Fauve J3| that by tuning the ratio of wavenumbers of the most unsta- 
ble modes one could control the resonant triad interactions that are important to the pattern 
formation process [^2| . Indeed all of the exotic patterns mentioned above were obtained with 
experimental parameters near the so-called "bicritical point" where two modes lose stability 



2 



simultaneously ||, ||, ||] . 

The role of resonant triad interactions in the formation of Faraday wave patterns has been 
investigated extensively by Vihals and co-workers for both the case of single frequency forc- 
ing |2^, [24|] , and two-frequency forcing 0] . The most detailed two- frequency calculations focused 
on the situation where the frequency ratio was 1/2, and the onset surface wave response was 
subharmonic with the forcing. Zhang and Vihals compared their theoretical results with exper- 
imental results of Muller Q, who observed subharmonic hexagons, triangles and squares near 
the bicritical point; which pattern was observed depended on a relative phase between the w 
and 2tt> sinusoidal waveforms in the forcing function. A key theoretical idea behind the work 
of Zhang and Vihals jj], Q , and Chen and Vihals 1 23 is that the presence of certain resonant 
triads composed of Fourier mode wavevectors ki, k 2 and ki + k 2 (|ki| = |k 2 |) can suppress 
the formation of regular wave patterns that involve the ki and k 2 modes. In particular, this is 
the case when the ki and k 2 modes are excited at the onset of the Faraday wave instability, 
while the mode with wave vector ki + k 2 is only weakly damped. This situation arises quite 
naturally near the bicritical point in parameter space. 

Motivated by the ubiquity of quasi-patterns in two-frequency Faraday experiments, Lifshitz 
and Petrich || recently investigated twelve-fold quasi-patterns within the framework of a simple 
Swift- Hohenberg type model for the evolution of a real scalar field u(x, t). Their model lacks 
reflection-symmetry u — > —u, and the linearized equations lead to a neutral stability curve 
with a double absolute minimum. However, because their model leads to time-independent 
patterns arising from steady state bifurcation, it does not capture one of the key features of the 
two-frequency Faraday problem. Specifically, at the bicritical point in the Faraday problem, 
one of the neutral curve minima is associated with a Floquet multiplier (FM) +1, while the 
other minimum is associated with a FM of —1. In this paper we are specifically interested in 
this feature of the bicritical point of the two- frequency Faraday experiment, and its potential 
influence on the resonant triad interactions important to the pattern formation problem. 

We investigate the role of resonant triads in the formation of wave patterns near the bi- 
critical point in parameter space. We focus on the situation that applies to the two-frequency 
quasipattern || and superlattice pattern g experiments in which the forcing frequency ratio, 
m/n < 1, has m even and n odd. We show that the usual contribution of resonant triads to 
the cross-coupling coefficient in the amplitude equations may be greatly suppressed in this case. 
This is in marked contrast to the situation of subharmonic waves with forcing ratio m/n < 1, 
where m is odd and n is even, e.g. m/n = 1/2, as studied by Zhang and Vihals. In this latter 
case, resonant triads lead to peaks in the cross-coupling coefficient g(9) (|]. In order to under- 
stand the distinction between these two situations, we first recall that in each of these cases 
the neutral stability curve has a double absolute minimum at the bicritical point, with one of 
the minima associated with a FM= +1 and the other with a FM= —1 [5L m]. At the bicritical 
point, there are two distinct critical wavenumbers, k m and k n (k m < k n ). The FM= +1 is 
associated with the smaller critical wavenumber k m if the even frequency is less than the odd 
one (i.e., m/n = even/odd); and it is associated with the larger wavenumber k n if the even 
frequency is greater than the odd one (i.e., m/n — odd/even). We show that in the former 
case the normal form of the coupled amplitude equations describing the relevant resonant triad 
interaction does not possess any quadratic terms because of a symmetry associated with the 
subharmonic waves. This is in contrast to the case of odd/even forcing where quadratic terms 
are present in the normal form. 
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In this paper, we demonstrate this distinction between even/odd and odd/even forcing 
by considering a simple example of one-dimensional waves that are parametrically excited by 
two frequencies. Rather than consider the full hydrodynamic equations, we investigate this 
issue using a simpler model derived from the free-surface Navier-Stokes equation by Zhang and 
Vinals pi]; their model applies to a deep layer of low- viscosity fluid. We consider both 2w/3w 
and \u /2lh forcing frequencies, over a range of frequencies u>, including a critical frequency for 
which a one-dimensional spatial resonance occurs. Specifically, at the critical frequency, the 
minima at the bicritical point are in ratio k m /k n = 1/2. We find that this spatial resonance 
leads to a divergence of the Landau coefficient in the amplitude equation describing bifurcation 
to subharmonic waves in the case of the 1/2 forcing frequency ratio; this is consistent with the 
results of Zhang and Vinals j24j. In contrast, the Landau coefficient for waves, parametrically 
excited by two-frequency forcing in ratio 2/3, is unaffected by any parameter proximity to the 
bicritical point and/or the spatial resonance point. While the wave vectors of the critical modes 
ki and 2ki are in resonance, there is a mismatch in their frequencies: in this case the usual 
large contribution to the Landau coefficient due to the spatial resonance is absent. 

The next section of the paper provides the necessary background to our analysis. It reviews 
the key theoretical ideas about the role of resonant triad interactions in pattern formation prob- 
lems, as well as the linear theory for Faraday waves, parametrically excited by two-frequency 
forcing. This background section closes with a discussion of the normal form symmetries as- 
sociated with the subharmonic instability. Section || sets up the example of one-dimensional 
surface waves, indicating the precise analogy between triad resonance in two dimensions and 
spatial resonance within this simpler one-dimensional problem. This section also presents the 
governing hydrodynamic equations that describe parametrically excited surface waves in the 
limit of an infinite depth fluid layer and weak viscosity. The last part of Section || sets up the 
weakly nonlinear analysis that we use to compute the cubic Landau coefficient for harmonic 
and subharmonic waves. Section ^| presents bifurcation results for two examples: waves para- 
metrically excited by two-frequency excitation in the ratio 2/3 and waves excited with a 1/2 
excitation frequency ratio. Section |g summarizes our results, and indicates some directions for 
future work. 



2 Background. 

2.1 Resonant triad interactions. 

One of the central ideas in pattern formation studies of quasi-patterns || |§|, |l4|, [2^], and 
surface wave patterns parametrically excited by two-frequency forcing 0, || is that of resonant 
triad interactions. In particular, there is interest in systems, such as the Faraday system with 
two-frequency forcing, for which it is possible to have a neutral stability curve with a double 
minimum. The resonant triad interaction of interest is dictated by the locations of the minima, 
say k m and k n , of the neutral curve. The resonant triad is made up of wave vectors ki, k2 and 
ki + k 2 , where |ki| = k 2 = k m and |ki + k 2 = k n . Thus the angle 6 between ki and k 2 
satisfies 




Here we have assumed that < k m < k n < 2k m ; see Figure [j]. 
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Figure 1: On the left is a plot of a neutral stability curve A(k), showing minima at k — k m and 
k = k n . On the right is an associated resonant triad ki, k2 and ki +k2, where |ki| = | k2 1 = k m 
and |ki + k 2 | = k n . The angle 9 between ki and k 2 is related to the ratio k n /k m by (|l|). 



To illustrate the possible influence of resonant triad interactions on pattern formation, we 
first consider a steady state bifurcation problem, where modes of wavenumber k m and k n 
lose stability (almost) simultaneously with the increase of an external control parameter. We 
consider a bifurcation problem on a six-dimensional center manifold for linear modes 

A(t)e lkl - X + B(t)e*= x + C(t)e l(kl+k2) - x + c.c, (2) 

where A,B,C £ C. Symmetry considerations determine that the unfolding of the bifurcation 
problem takes the form, through cubic order in the amplitudes, 

A = XA + aB*C+(a\A\ 2 +b\B\ 2 + c\C\ 2 )A 

B = \B + aA*C+ (a\B\ 2 + b\A\ 2 +c\C\ 2 )B (3) 
C = [iC + f3AB + (d\A\ 2 + d\B\ 2 + e\C\ 2 )C , 

where the asterisk indicates the complex conjugate, and the coefficients are all real. If A = 0, 
/i 7^ 0, then one can further reduce the bifurcation problem to one involving the critical modes 
A and B, with C constrained to the center manifold: C = —-AB + • • •. One obtains, for |A| 
sufficiently small, the reduced (unfolded) bifurcation problem 

A = XA + a\A\ 2 A+ (b- ^j\B\ 2 A 

B = XB + a\B\ 2 B + (b-^\A\ 2 B . (4) 

The presence of the near critical mode C in Q leads to a large cross-coupling coefficient 
g(6) = (b — 2^) in the amplitude equations (Q) since \/j,\ C 1. A consequence of this is that 
patterns that involve an equal amplitude superposition of modes A and B tend to be unstable 
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at onset. For example, the stability of steady rhombic patterns at angle 8 (A = B in (0)), 
within the setting of the amplitude equations ([I]), is determined by two eigenvalues whose signs 
are sgn(a + b — ^) and sgn(a — b + ^). When |^| ^S> \a\ + \b\, the two eigenvalues take 
on opposite signs and the pattern is necessarily unstable at onset. If the spatial resonance 
occurs near 9 = tt/2, for example, then square planforms are unstable. Quasi-patterns, and 
spatially-periodic "superlattice" patterns, such as described in J|, [| involving both of the 
modes A and B above are similarly destabilized by a large cubic cross-coupling coefficient g(8) 
in the amplitude equations. This point is discussed for various regular patterns, including 
quasi-patterns, in 0]. 



2.2 Bicritical curves for the two-frequency Faraday problem. 

We now examine, in greater detail, the double minimum of the neutral curve associated with the 
Faraday problem with two-frequency forcing. Besson, Edwards and Tuckerman Q computed 
the linear stability of the flat free surface of a shallow fluid layer subjected to a periodic vertical 
acceleration 

g(t) — go + <7 z [cos(x) cos(mujt) + sin(x) cos(nujt + </>)]. (5) 

Here go is the usual gravitational acceleration, m and n are co-prime integers, and g z , u>, x, <t> are 
additional external control parameters. Note that the angle x controls the relative amplitudes 
of the mw-forcing and nw-forcing. Besson et al. [0 compared, with good agreement, their 
theoretical predictions with their experimental results for various fluids with different viscosities 
in the case that m = 4 and n = 5 in §). In this case, they found that the initial instability 
of the fluid surface was either associated with a Floquet multiplier (FM) —1 or +1, depending 
on whether the odd or even frequency component in (|^) dominated. The transition between 
subharmonic (FM= —1) and harmonic response (FM= +1) occurred at a particular value of 
the angle x m ©i called the bicritical point. 

Here we do not consider the full hydrodynamic problem, but rather a simplified model 
derived by Zhang and Vinals [Q, which applies to a deep layer of weakly viscous fluid. The 
linearized equation for the free surface mode hk(t)e lkx is 



til + Avk 2 ti k 



g(t)k + + W h k = 0, (6) 

P J 



where v is the kinematic viscosity of the fluid, T is the surface tension, and p is the fluid density. 
Note that in the absence of the time periodic parametric forcing (i.e. g z = 0), one recovers 
a dispersion relation, where the damping is given by 2vk 2 and the water wave frequency u>o 
satisfies the usual gravity-capillary wave dispersion relation Uq = g$k + Tk 3 / p. 

In Figure ^ we present an example of the neutral stability curves g z [k) for the parameters 
used by Besson, et al. Q, and for various values of the parameter x- These curves were obtained 
from (^|) using the same method that Besson, et al. Q used to obtain the neutral curves for 
the full hydrodynamic problem. The bicritical point in this example occurs at x — 63.11°, for 
g z = 3.13 go (go = 980.665 cm/s 2 ), k + = 7.14 cmT 1 and fc_ = 9.35 cm -1 , where k± are the 
values of k associated with FMs ±1. The locations of the minima can be roughly approximated 
by considering the simple water wave dispersion relation uj 2 = gok + Tk 3 / p, and assuming 
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Figure 2: Neutral stability curves computed from (|6J). The vertical axis is g z /go and the 
horizontal axis is k; the angle \ m (H) 1S indicated in the lower right corner of each plot. Floquet 
multipliers of +1 (—1) arc indicated by solid (dashed) lines. The first instability encountered 
with increased forcing g z is to harmonic waves if x < Xc ~ 63°, and to subharmonic waves 
for x > Xc- The other parameters of the forcing function are m — A, n = 5, 4> — and 
u}/2tt = 11 Hz; the fluid parameter are T — 20.6 dyn/cm, v — 0.209 err? j s and p = 0.95 g/cm 3 . 
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subharmonic response to the forcing frequencies at mui and ntv, i.e., for m = 4 and n = 5: 



4w\2 . rfci 



5o 



-5w\2 , rfc 3 



(t) m **- + — (7) 

For the example of Figure [2] this approximation yields k + w 8.0 and fc_ w 9.8. We note that 
the neutral stability curves obtained from the simple linear problem (^) are qualitatively very 
similar to those obtained by Besson, et al. from the full hydrodynamic problem, cf. our 
Figure |^ with Figure 2 in j^j . 

2.3 Normal form symmetries. 

We now re-examine the role of resonant triad interactions for the mode interaction problem 
pertinent to the two-frequency Faraday experiment. In particular, we will investigate the 
additional effect of normal form symmetries on the pattern formation problem. We wish to 
contrast the situations where the forcing frequencies muj and nco in (^), m < n, have m even, 
n odd with m odd, n even. We are specifically interested in the case where the angle x hi (|^) 
is close to the bicritical point Xc- 

We analyze the resonant triad interaction in terms of a stroboscopic map since we are 
interested in a periodically- forced system. We denote the surface height at time t — pT, p an 
integer, T = 2tt/lu, by h p (x). Let the resonant triad be 

h p = Ape 4311 " + B p e lk2X + C p e i( - kl+k2 >* + c.c. + • • • , (8) 

where |kx| = |k 2 | = k m and |kx + k 2 | = k n , < k m < k n < 2k m . Here the m,n subscripts 
indicate that the critical wavenumbers can be roughly associated with the muj and ruv forc- 
ing, respectively, as in (0) above. The general form of the cubic-order amplitude equations, 
consistent with translation symmetry and reflection symmetry, is (cf. equation ||) 

A p+1 = XA p + aB* p C p + (a\A p \ 2 + b\B p \ 2 + c\C p \ 2 )A p 

B p+1 = XB P + aA* p C p + (a\B p \ 2 + b\A p \ 2 + c\C p \ 2 )B p (9) 

C p+1 = iiC p + I3A P B P + (d\A p \ 2 + d\B p \ 2 + e\C p \ 2 )C p . 

In the case that m is even and n is odd, the Floquet multipliers for the linearized problem are 
A = +1 and fi = — 1 when x — Xc, whereas if m is odd and n is even, then A=— 1, /i = +l. 

It may be possible to further simplify the bifurcation problem (||) by a normal form trans- 
formation. In particular, there exists a near-identity nonlinear transformation such that all 
nonlinear terms in (|9|), which do not commute with the matrix L of the linearized problem, can 
be removed (see, for example, Crawford pq|). Here 



(10) 
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where |A| = \fi\ = 1. 

In the case that A = — 1, /i = +1, i.e., m odd and n even in the forcing ([&]), the bifurcation 
problem (^|) is already in normal form since the normal form symmetry is equivalent to a 
translation by d, where ki ■ d = k2 ■ d = 7T. The presence of the quadratic terms in this normal 
form means that the resonant triad interactions will have a strong influence on the pat tern 
formation problem, as described for the simple steady state bifurcation example in Section 2.1. 

In contrast, in the case that A = +1, fi = — 1, i.e., m even and n odd, the normal form 
transformation allows the quadratic terms in the bifurcation problem to be removed. The 
normal form of the bifurcation problem, through cubic order, is 

A p+1 = A p + {a\A p \ 2 + b\B p \ 2 + c\C p \ 2 )A p 

Bp+i = Bp + (a\B p \ 2 + b\A p \ 2 +c\Cp\ 2 )B p (11) 
C p+1 = -C p + (d\A p \ 2 + d\B p \ 2 + e\Cp\ 2 )C p . 

In this case, the bifurcation to harmonic waves is investigated in the invariant subspace C = 0, 
and there is no divergence of the cross-coupling term due to the resonant triad interaction. 
The usual contribution of the resonant triad interaction to the pattern formation problem 
is suppressed in this case. In the next section we test these ideas by performing explicit 
computations of bifurcation coefficients from the simplified hydrodynamic model, due to Zhang 
and Vihals [E3, of the two- frequency Faraday experiment. 



3 Two-frequency forcing of one-dimensional surface waves. 

3.1 Resonant interactions in a one-dimensional problem. 

It is possible to investigate the effects of strong spatial resonance on a pattern formation problem 
in one spatial dimension. This is done by considering the situation where the minima k m and 
k n of the neutral curve satisfy k m /k n = 1/2. In this case, the stroboscopic map associated with 
the critical modes Z p e lkx + W p e 2lkx + c.c. is, through cubic order, 

Z p+1 = XZp + aWpZ; + (a\Zp\ 2 +b\Wp\ 2 )Z p 

W p+1 = fiWp + (3Z 2 p + (c\W p \ 2 + d\Z p \ 2 )Wp . (12) 

In form, this is identical to (^) restricted to the subspace A = B = Z, C = W. 
If A = —1 and fj, = +1, then ( |l2| ) is already in normal form. The unfolding is 

Zp+i = -(l + e)Zp + aWpZ; + (a\Z p \ 2 + b\Wp\ 2 )Z p 
W p+1 - {l + S)Wp+pZ 2 + {c\W p \ 2 +d\Z p \ 2 )Wp . (13) 

Subharmonic waves bifurcate at e = 0. For S ^ 0, the center manifold is given by W = 
—f3Z 2 /8 + • • •, and the branching equation is then given by 

= -e+{a~a[3/5)\Z\ 2 . (14) 

We see that the presence of the near critical spatially-resonant harmonic mode W forces the 
subharmonic wave to bifurcate with very small amplitude, since the cubic coefficient is large 



9 



when 5 is small. This effect is the one-dimensional analogue of the cross-coupling term g(ff) 
diverging near a spatial resonance. This is the situation that can occur when m is odd and n 
is even (n > m) in (^), e.g. m = 1, n = 2. 

If, on the other hand, A = +1, /i = — 1, then the quadratic terms in ( |T^ ) can be removed, 
and the usual effect of the spatial resonance on the pattern formation problem is suppressed. 
Specifically, the (unfolded) normal form of the bifurcation problem is 

Z p +i = {l + e)Z p + (a\Z p \ 2 +b\W p \ 2 )Z p 

W p+ i = -{l + S)W p + {c\W p \ 2 +d\Z p \ 2 )W p . (15) 

Harmonic waves, which bifurcate from the trivial solution at e = 0, are investigated in the 
invariant subspace W = 0; they satisfy the branching equation 

= e + a\Z\ 2 . (16) 

Subharmonic waves bifurcate from the trivial solution at 8 — and can be investigated in the 
invariant subspace Z = 0; they satisfy the branching equation 

= -5 + c\W\ 2 . (17) 

Neither of these branches is affected by the spatial resonance. This is the situation we expect 
when m is even and n is odd in (§), e.g. to = 2, n = 3. 

In Figure |^, we plot the bicritical point Xc as a function of lo for both (m, n) = (2, 3) 
and (m,n) = (1,2). These curves apply to the linear equation (||), with T = 20.6 dyn/cm, 
v = 0.209 cm 2 1 s and p = 0.95 gjcrr? . In this figure we show also the ratio of wavenumbers 
kn/km at the bicritical point. Note that in each case the spatial resonance k n /k m — 2 is 
achieved for a particular value of the forcing frequency lo — uj res . 

3.2 The Zhang- Vinals hydrodynamic model. 

Zhang and Vinals Q derive the following model for the surface deviation h(x, t) and surface 
velocity potential $(x, t), where x 6 R 2 : 

{d t - 7V 2 )/i - 2>$ = T(h, $) (18) 
(d t ~ 7V 2 )$ - Ah = G(h, $) . 

Here the linear operator A is 

A = f o V 2 - G - 4/[sin(2i) + k sin(2pt + 0)] , (19) 

where p = n/m. The linear, nonlocal operator T> multiplies each Fourier component of the field 
by its wavenumber modulus. T and Q are nonlinear operators to be defined shortly, together 
with the parameters 7, r , Go, / and k. 

In deriving the above, Zhang and Vinals assumed that the experimental forcing takes the 
form <7 z r(sin(2mcL;o£) + K sin(27iijJo^ + </>)), where k = (1 — r)/r; they have scaled time by 
the frequency tjiloq. Here we choose a slightly different scaling of time; we assume that the 
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bicritical line, m=1,n=2 bicritical line, m=2,n=3 




70 1 1 ' 1 1 1 1 1 1 1 1 1 1 1 1 
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Figure 3: Bicritical point \c plotted as a function of the frequency w in @, where r/> = and 
r = 20.6 dyn/cm, v — 0.209 err? j s and p = 0.95 g/cm 3 in (||). Plotted below each bicritical 
line is the ratio of critical wavenumbers k n /k m for x = Xc(^). On the left m = 1, n = 2, and 
on the right m = 2, n = 3. 
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parametric forcing takes the (dimensionless) form /[cos(x) cos(tot) + sin(x) cos(nr + <f>)]. We 
let t = mr/2 and $ = — $ in (ITsh to obtain 



(d T - 7 V 2 )/i - X>$ = F(h,<f>) (20) 
(0r - 7V 2 )$ - .4A = 

where 

7-?^. (2D 

The wavenumber &o satisfies the equation 

9ok + ^ = ——, 22 
p 4 

where w = 2wo is the dimensioned forcing frequency. The linear operator A is now 

A = T V 2 - G - /(cos(x) cos(mr) + sin(x) cos(nr + </>)), (23) 

where 

r _ rfc o r _ ff fc o . ff z fco , 9 .s 

The nonlinear operators are 

jf(/i,$) = -v-(/iV$) + iv 2 (/i 2 i5$)-p(/ii5$) + 25[/iP(/iP$) + i/i 2 v 2 $] (25) 
0(/i,$) = i(5$) 2 - i(v$) 2 - (p$)[w 2 $ + p(/ip$)] - ir v • (v/i(v/i) 2 )). 

Finally, we put the governing equations in a convenient form for the perturbation analysis 
that follows. In order to recover the linearized equation (||) for h, we apply the operator 
(d T — 7V 2 ) to the first of equations (|2p|), and then use the second equation to express (9 T — 7 V 2 )<I> 
as Ah + Q(h, 5>). We supplement this equation with the first of equations (po|), re-written as 
an equation for the field to obtain the full system of equations: 

Ch = (d r -jV 2 )T(h,<S>)+Vg(h,<i>) (26) 
X>$ = (d r - 7 V 2 )/i- H h ,®) ■ 

Here the linear operator is 

C = d TT - 2 7 V 2 <9 r + ( 7 2 V 4 - VA) , (27) 
where the operator A is given by (p3|). 

3.3 Weakly nonlinear analysis: derivation of the bifurcation problem. 

We restrict our analysis to spatially one-dimensional solutions of the variable ieR, and use a 
two-timing perturbation method to determine weakly nonlinear solutions of the system (pq) in 
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the vicinity of the bifurcation to (sub)harmonic waves, i.e. near a primary instability associated 
with a Floquet multiplier +1 (—1). To do this we introduce a small parameter e, such that 



h(x,r) = eh 1 (x,T,T) + e 2 h 2 (x,T,T) + e s h 3 (x,T,T)+--- (28) 
$(ar,T) = e$ 1 {x,T,T) + e 2 $ 2 (x,T,T)+e 3 $ 3 (x,T,T) + --- , 

where 

T = e 2 r , f = f + e 2 f 2 . (29) 

Here /o is the value of the forcing at the bifurcation point. We seek spatially-periodic solutions 
in the following separable Floquet-Fourier form: 

h x = Zl {T) Pl {T)e^ x + c.c. h 2 = zl{T)p 2 {T)e i2k " x + c.c. (30) 

$i = Zl (T) qi {T)e^ x + c.c. $ 2 = zl(T)q 2 { T y 2k ° x + c.c. 

Here pi, qi (I — 1, 2) are 27r-periodic functions of the fast time r in the case of harmonic waves 
(FM=+1); in the case of subharmonic waves they are 47r-periodic in r. The wavenumber k c is 
associated with the first unstable mode. 
At leading order in e, we obtain 

C hi = (31) 
P$! = {d T -~td xx )hi 

Here the linear operator Cq is C with / = /o, where C is defined by (p7|). The equation £o^i = 
determines the bifurcation point /o by the solvability condition that it have a periodic solution; 



it is the equation we solved in Section 2.2 to determine the neutral stability curves. This 
equation also determines the critical wavenumber k c by the requirement that /o be the smallest 
value of / that admits a periodic solution. Given /o and fc c , we determine the solution of 
Cohi =0 involving the periodic function £>i(t), which we may assume is real. We express p\ 
in terms of its (truncated) Fourier series 

N 

Pi W =^a j e i ( j+ ^ T + c.c, (32) 

3=0 

where N is chosen large enough so that the solution is well converged. Here [i = for the 
case of harmonic waves and /i = 1/2 for subharmonic waves. We readily solve the second of 
equations ( |3l] ) for $i, yielding 

q 1 (r) = ^(d T + 1 k 2 c )p 1 . (33) 

At order e 2 , we obtain 

C h 2 - -(a T - 7 9 M )[&(/ii&$i)+P(/ l iP$i)] + ip[(P$i) 2 -(&$i) 2 ] (34) 



2 



2fc3 z 2 g 2 e 2lfc ^ + c.c. 
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where q\ is a real 27r-periodic function. We determine a solution h 2 involving the real periodic 
function p 2 (r), given by its (truncated) Fourier series representation 



N 

P2 = ^bje 1 ^ + c.c. 
j=o 

We also find, at this order, that $2 satisfies 

v$ 2 = (d T -^d xx )h 2 

since d x (hid x $i) + V{hiV<$>i) = 0; thus 



12 (r) = — (d r + 4 7 fc c 2 )p 2 . 



(35) 



(36) 



(37) 



Finally, we consider the equation at order e 3 : 

£ ^3 = — 2d T dxhi + 2^d xx dTh\ — /2[cos(x) cos(mT) + sin(x) cos(nr + (f>)]Vhi 

+ (d T - ld xx ) \-d x {hid x <S> 2 + - V{hiV<S> 2 + h 2 V<S>i) (38) 



-\d xx {h\v$i) + v[h 1 v(h 1 v<$> 1 ) + ^h\d xx ^] 



V 



(D$i)(D$2 



-(9 x $i)(a x $ 2 ) - S*iMx**i + £(*n£*i)] - -y9x(9x/n) 3 

In order to ensure that a 27r-periodic solution exists, we must apply a solvability condition to 
this equation, written compactly as £0^3 = H(x,t,T). Specifically, we require 



(hi,£ h 3 ) = (h\, H(x,t,T)) = 0, 



where the inner product is 



hi = pi(r)e cX is a periodic solution to the adjoint linear problem C' Q hi = 0. Here 

C\ = (d T + -fd xx ) 2 - VA, 

where / = /o in A. 

The solvability condition leads to the amplitude equation 

S^=af 2 zi + (A + B)\z 1 \ 2 zi, 



(39) 



(40) 



(41) 



(42) 



where 



1 f i7T 

- tJ « 



+ z^ipi + jk 2 pi)pi dr 
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a = — - / [cos(x) cos(mr) + sin(x) cosfnr + </>)]piPi dr (43) 

A = -7T {{qm)' + ifJ.qxP2 + jk^qiP2)px dr 
27r Jo 

k 3 f iw 3 
B = ^-J (-(pfg-i)' - W?<7i ~ lk 2 c p\qi + k c q\p x + 2 k c T o\pi\ 2 Pi)pi dr. 

In the above, differentiation with respect to r is denoted by a prime. We also took the periodic 
solution pi to be real. 

In the above derivation of the bifurcation problem (|4^), we have separated the two contri- 
butions A and B to the cubic coefficient A + B. The contribution A comes from the quadratic 
nonlinear terms in the original hydrodynamic model; it depends on the spatially resonant modes 
e i2k c x ag j g gyjdent from the P2 terms in the integral expression for A given in (^3|). The con- 
tribution B comes from cubic nonlinearities and therefore depends only on the mode e ikcX . We 
calculate all of the coefficients 6, a, A and B numerically as follows. 

At leading order in e, equation ([n]) for hi reduces to the linear problem (Q). The solution 
for hi is found by representing Pi(t) as the truncated Fourier series given in equation (|32[). As 
discussed in , the problem then reduces to a generalised eigenvalue problem of the form 

Ax = g 2 Bx, (44) 

where x is a vector of the Fourier coefficients a,j in ( |32| ) . The EISPACK routine rgg is used to 
find the eigenvalues and corresponding eigenvectors, the eigenvalues giving the linear stability 
curves shown in Figure 0. The eigenvalue corresponding to the minimum of the resonance 
tongue under consideration is then determined, with the corresponding eigenvector then giving 
the coefficients a,j for the spectral representation of Pi(t), from which qi(r) can be calculated. 
The adjoint linear problem for pi is solved in a similar way to that for p\. 

Once a spectral representation for qi is found, q\ is calculated using a pseudospectral ap- 
proach. This then enables the order e 2 problem for /12, equation (|34|), to be written as a 
nonhomogeneous linear problem for pi- This can be solved using the EISPACK routine rg. Fi- 
nally, the coefficients 6, a, A and B in the amplitude equation (^2|) are computed by calculating 
the various products of pi , P2 and qi using a pseudospectral method and then calculating the 
inner products. 

Typically 20 Fourier modes sufficed to represent pi,P2 and qi, and 257 collocation points 
were adequate in the pseudospectral calculation of the nonlinear terms. Checks were done, with 
twice as many modes and collocation points, to ensure that the results were well converged. 

In the next section we present plots of the ratio A/B of the two contributions to the cubic 
Landau coefficient for the case of bifurcation to waves excited with forcing frequency ratio 2/3, 
and for the case of bifurcation to waves excited with forcing frequency ratio 1/2. 

4 Results. 

We have calculated the coefficients 5, a, A and B in the amplitude equation (^) for four cases. 
These correspond to the two lowest resonance tongues in the two cases m = 1, n = 2 and 
m = 2, n = 3: whether harmonic or subharmonic waves bifurcate at a smaller value of the 
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m=2,n=3 harmonic waves 



m=2,n=3 subharmonic waves 




10 20 30 40 50 60 70 80 90 100 10 20 30 40 50 60 70 80 90 100 

omega omega 

Figure 4: Plots of A/B as a function of \ an d u>, for m = 2, n = 3 and = in (||). 
Computations are done at the bifurcation point to (sub)harmonic waves on the left (right). 
The fluid parameters are T = 20.6 dyn/cm, v = 0.209 cm 2 / s and p = 0.95 g/cm 3 . We give 
both surface plots (top) and the corresponding contour plots (below). In the contour plots we 
reproduce, from Figure the bicritical line in the w, x~plane. The square on this curve gives 
the value of oj where the 1/2 spatial resonance occurs. 

driving amplitude depends on whether uj and x take values which are above or below the 
bicritical lines given in Figure [| The other parameters used are the same as those listed in the 
caption to Figure |[ 

Figure^ shows the ratio A/B of the two contributions to the cubic Landau coefficient for the 
excitation frequency ratio 2 /3 for the two tongues which bifurcate at lowest amplitudes of the 
parametric excitation, g z . In each case the upper graphs show the surface representing A/B as 
a function of to and x- The lower graphs show the corresponding contour plots. Superimposed 
on the contour plot is the bicritical line from Figure |3| above: the square mark on this graph 
indicates the point on the bicritical line at which the ratio between the wavenumbers for the two 
minima is 2. Note that in this case, where the excitation frequency ratio is 2/3, the bifurcation 
to the harmonic waves occurs at lower amplitudes of the parametric excitation g z for values of 
oj and x below the bicritical line, while the bifurcation to subharmonic waves occurs first for 
values of ui and x above the bicritical line. The most striking feature of the plots associated 
with the harmonic waves in Figure |] is the insensitivity of the quantity A/B to the point on 
the bicritical line at which the ratio between the wavenumbers is 2. Spatial resonance is also 
not important for the subharmonic waves since the subharmonic tongue occurs at the higher 
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m=2,n=3 harmonic waves 




200 
omega 



200 
omega 



Figure 5: Plots of A/B similar to those in Figure [|, but computed for m = 1 and n = 2. The 
divergence of A/B for subharmonic waves is discussed in the text. This divergence occurs on a 
curve that is tangent to the bicritical line at the point ui res — 191 s _1 , Xc = 72.83°; this is the 
point where there is a 1/2 spatial resonance on the bicritical line. 
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Neutral stability curves, 111=1,11=2, omega=135, chi=77.8 
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Figure 6: Neutral stability curves for to — 135 s _1 ,x — 77.8°. The harmonic (subharmonic) 
tongue is indicated by a solid (dashed) line. The minimum of the subharmonic tongue occurs 
at k m = 3.308 cm -1 , when g z — 2.1658 go. Also marked is the point where 2k m intersects the 
harmonic tongue. Note that this intersection occurs at the critical value of g z , where bifurcation 
to subharmonic waves occurs. 



wavenumber in this case of even/odd excitation frequencies. This is in marked contrast to the 
case of odd/even excitation, as demonstrated by the graphs shown in Figure [5], which were 
computed for an excitation frequency ratio 1/2. This time, it is the subharmonic instability 
which occurs first for values of u> and x below the bicritical line and the harmonic instability 
which sets in first for values of u> and x above the bicritical line. As for the previous graph, 
A/B for the harmonic waves is unaffected by the ratio between the wavenumbers and the 
parameter proximity to the bicritical line. However, in the case of subharmonic waves there is 
a singularity occurring on the bicri tical line at the point of spatial resonance. This is consistent 
with the discussion in Section 



3.1 



above, showing how critical the nature of the instability is 
in determining whether or not spatial resonance is important to the pattern formation process 
near the bicritical line. 

Note that in Figure || there is in fact a line of singularities tangent to the bicritical line at 
the point on the bicritical line where the ratio between the wavenumbers is 1/2. This line of 
singularities can also be understood through spatial resonance. Recall that above the bicritical 
line the subharmonic instability occurs at a higher value of the excitation amplitude than the 
harmonic instability, as shown in Figure || for the values u> — 135 s" 1 ^ — 77.8°. These 
values were chosen as representative of a point on the line of singularities shown in Figure |^. 
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In this case the minimum of the subharmonic tongue occurs at an excitation amplitude of 
g z = 2.1658 go and at a wavenumber of k m — 3.308 cm" 1 . However, at this excitation amplitude 
there is a bifurcation to harmonic waves occurring at the wavenumber of 6.616 cm -1 = 2k m , 
and therefore spatial resonance will occur. Similarly, for values of u> and x to the right of 
the singularity which occurs on the bicritical line, spatial resonance can occur between the 
minimum of the subharmonic curve and the right side of the harmonic tongue. Note that, in 
practice, it is only close to the bicritical line that one would expect to observe the effect of the 
spatial resonance on the subharmonic waves, as in all other cases bifurcation to harmonic waves 
occurs at significantly lower amplitudes of the excitation. We also note that for parameters on 
the bicritical line the center manifold reduction described in Section 3.1 can break down. 



5 Conclusions. 

In this paper we have investigated how normal form symmetries affect the role of resonant triads 
in the pattern formation problem for surface waves parametrically excited by two-frequency 
forcing. We focused on the behavior of the system near the bicritical point in parameter space, 
where modes of wavenumber k m and k n lose stability simultaneously with one mode associated 
with a Floquet multiplier of +1 and the other associated with a Floquet multiplier of —1. Our 
analysis shows that when the Floquet multiplier +1 is associated with the smaller wavenumber, 
then quadratic terms may be removed from the relevant amplitude equations. In this case, the 
contribution of resonant triads to the bifurcation problem is not affected by proximity to the 
bicritical point. In contrast, when the Floquet multiplier —1 is associated with the smaller 
wavenumber, then the quadratic terms cannot be removed by a normal form transformation. 
Hence, in this instance, resonant triads influence greatly the bifurcation problem near the 
bicritical point. The former situation applies when the forcing frequency ratio is m/n < 1, 
with m even and n odd, while the latter situation occurs when m is odd and n is even. Thus we 
expect normal form symmetries are important to understanding the experimentally observed 
quasipatterns and superlattice patterns || which employ even/odd forcing. Such effects are 
necessarily neglected in theoretical models, such as the model of Lifshitz and Petrich in 
which quasipatterns form through a steady state bifurcation. 

This paper demonstrates the influence of the normal form symmetries on the bifurcation 
problem by considering the example of one-dimensional waves, parametrically excited by two- 
frequency forcing. In this one-dimensional problem a spatial resonance, involving Fourier modes 
of wavenumber k m and 2k m , takes the place of a resonant triad. Specifically, we focused on 
the situation where the bicritical point involves modes k m and k n = 2fc m . We considered 
instability to harmonic waves and subharmonic waves when the two-frequency forcing was in 
ratio m/n — 1/2 and m/n = 2/3. Rather than perform the weakly nonlinear analysis on the 
full hydrodynamic equations, we used the simpler Zhang and Vinals model [|l] that applies to a 
deep layer of a nearly inviscid fluid. Consistent with our general bifurcation theoretic analysis, 
we found that only in the case of subharmonic waves, parametrically excited by m/n — 1/2- 
forcing, did the presence of the bicritical point lead to a diverging Landau coefficient in the 
bifurcation problem. In the other cases, the Landau coefficient was completely insensitive to 
any parameter proximity to the bicritical point or to the spatial resonance. 

Our work suggests that the spatially resonant triads, important to a broad class of pattern 



19 



formation problems, do not play an important role in pattern formation of parametrically- 
excited surface waves near the bicritical point in the case of even/ odd forcing. However, since 
our analysis focused only on the case of a one-dimensional pattern formation problem, more 
analysis in the two-dimensional case is warranted. In the future we hope to carry out a more 
extensive analysis of the contribution of resonant triads to two-dimensional Faraday waves 
with even/ odd forcing. For instance, it would be of interest to carry out such an analysis for 
the physical parameters of the recent Faraday experiments in which superlattice patterns were 
observed ||. In these experiments the forcing frequencies were in ratio 6/7; the spatial Fourier 
transform of the harmonic wave patterns exhibited peaks, some of which could be associated 
with the two different frequency components of the forcing wave form ]5| . 

Although the analysis presented in this paper was motivated by and applied to the problem 
of Faraday waves with two-frequency forcing, we expect many of the symmetry-based ideas 
to carry over to other parametrically-excited pattern forming systems. We have in mind, for 
example, the recent experiments on one-dimensional surface waves on ferrofluids, which are 
excited by a time-periodic magnetic field |27j . In this system, both harmonic and subharmonic 
response occur with single-frequency forcing. 
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